arXiv: 1503.02147vl [math.CA] 7 Mar 2015 


Pade interpolation and hypergeometric series 


Masatoshi Noumi 

Department of Mathematics, Kobe University 


Abstract: We propose a class of Pade interpolation problems whose general solution 
is expressible in terms of determinants of hypergeometric series. 

Key words : Pade interpolation; hypergeometric series; Dodgson condensation; Krat- 
tenthaler determinant 

2010 Mathematics Subject Classification : 41A05; 33C20, 33E20 

1 Introduction 

In this paper we investigate a class of Pade interpolation problems to which the solutions are 
expressible in terms of determinants of hyper geometric series. Pade interpolation problems have 
been discussed by Spiridonov-Zhedanov p3] from the viewpoint of biorthogonal rational functions. 
They are also sources of the Lax pairs for discrete Painleve equations constructed by Yamada ESI, 
E3, and by NoumL-Tsujimoto-Yamada ED- The goal of this paper is to clarify how hyper geometric 
series arise in Pade interpolation problems, by analyzing the determinantal expression of the general 
solution. 

In Section 2, we formulate a general Pade interpolation problem and a universal determinant 
formula for the general solution (Theorem 12.IP We also show that the determinants expressing 
the general solution can be condensed to smaller determinants by a variation of the Dodgson con¬ 
densation (Theorem 12.2j) . After these preliminaries, we investigate in Section 3 a class of Pade 
interpolation problems relevant to generalized hyper geometric series r +i F r . We propose there two 
types of formulas expressing the solutions in terms of determinants of generalized hypergeometric 
series. The first one (Theorem 13. ip . derived through Theorem 12.21 is based on the condensation 
of determinants and Krattenthaler’s determinant formula, while the second (Theorem 13.21) is con¬ 
structed by means of the Saalschiitz summation formula for terminating 3 T 2 series. We remark 
that Pade approximations to generalized hypergeometric functions have been discussed by Luke 
E3, El- It would be an important question to clarify the relationship between interpolations and 
approximations in the context of generalized hyper geometric functions. Section 4 is devoted to the 
extension of these results to three types of very well-poised hypergeometric series including basic 
(trigonometric) and elliptic hyper geometric series. The two determinant formulas of Theorem 14.11 
and Theorem 14.21 are obtained by Warnaar’s elliptic extension of the Krattenthaler determinant 
and by the Frenkel-Turaev summation formula for terminating 10 V 9 series. 
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Two fundamental tools of our approach are the condensation of determinants along a moving 
core (an identity of Sylvester type), and variations of Krattenthaler’s determinant formula. For 
the sake of convenience, these subjects are discussed separately in Appendix A and Appendix B 
respectively. Generalization of Sylvester’s identity on determinants has been developed extensively 
by Miihlbach-Gasca [9] (see also El)- The version we use in this paper (Lemma IA.2I) . based 
on the Neville elimination strategy , is originally due to Gasca-Lopez-Carmona-Ramirez [2]. We 
also remark that Sylvester’s identity and its extensions play important roles in recent studies of 
integrable systems (see Spicer-Nijhoff-van der Kamp m for example). As to Appendix B, basic 
references are the works of Krattenthaler and Warnaar [15] (see also Normand [10] for recent 

works on the evaluation of determinants involving shifted factorials). Although the contents of 
these appendices are basically found in the literature, we include them as self-contained expositions 
which might be helpful to the reader. 

Throughout this paper we use the following notation of submatrices and minor determinants. 
For an m x n matrix X = { x ij)i <i<m 1<J<n (with entries in a commutative ring), we denote by 


Y"^l _ 


x hji • • • x hjs 


- X irjl ‘ ‘ ‘ X irjs - 



l<a<r,l<b<s 


( 1 . 1 ) 


the r x s submatrix of X with row indices i\,... ,i r £ {1,..., m} and column indices ji..... j r £ 
{1,... ,??,}. When r = s, we denote by det the corresponding minor determinant. 


2 Pade interpolation problems and their determinant solutions 

In this section we formulate general Pade interpolation problems and propose some universal de¬ 
terminant formulas for the solutions. 

Let /o(x), ..., fm(x) and go(x), gi(x), ..., g n (x) be two sequences of linearly independent 

meromorphic functions in x £ C and set N = m + n. We consider a pair (P m (x),Q n (x)) of two 
functions 


Pm(x') — PmO fo(x) T Pml f\{x) T ' ' ' T Pmm fm(.X ), 

Qn(x) = q n ogo(x ) + q n i9i(x) + ••• + qnng n (x), (2.1) 


which are expressed as C-linear combinations of fj(x) and gj(x) respectively. Noting that the ratio 
Pm(x)/Q n (x) contains N + l= m + n + \ arbitrary constants, we investigate the interpolation 
problem 


Pm{u o) P m («l) Pm{u]\f) 

Qn{u 0 ) V °’ Q n {ui) Vl ’ Q n {u N ) 


( 2 . 2 ) 


2 









for a set of N + 1 generic reference points x = uq, u\, ..., un and a set of N + 1 prescribed values 
vo, vi, ..., vn- This problem is equivalently rewritten as 

Pm{uk) ■ Qn(uk ) = A k ■ Pk {k = 0,1,... , N) (2.3) 

for v k = X k /p k (k = 0,1,, N). We remark that the Pade interpolation problem defined as above 
contains the Lagrange interpolation problem as a special case where n = 0 and go(x) = 1. 

A general solution of this Pade interpolation problem is given as follows in terms of (N + 2) x 
(N + 2) determinants: 


fo{x) 

fm(x) 0 

0 

Mo/o(«o) ••• 

hofm(u 0 ) X 0 go(u 0 ) ••• 

Xog n {u 0 ) 

APv/o(«iv) 

PNfm{uN ) XnPo{un ) 

XN9n{uN). 

0 

0 g 0 (x) 

9n{x) 

hofo{uo) 

Tofm(uo) Xogo(uo) 

Xo9n(u 0 ) 

-APv/o(“iv) 

•• PNfmiuN ) Xn9q{un) 

'• X N g n (u N ) 


(2.4) 


(2.5) 


Theorem 2.1 The pair of functions (P m (x),Q n (x)) defined by (12.41) . (12.51) solves the Pade inter¬ 
polation problem (E3D */ ( Pm(uk ), Qn(u k )) / (0, 0) for k = 0 , 1 ,... , N. 


In order to prove that this pair (P m (x),Q n (x)) gives a solution of the interpolation problem, 
we introduce two parameters A, p and consider the (N + 2) x (N + 2) determinant 


Rm, n (x;X,p) = det 


hfo(x) ••• hfmfx) Xg 0 (x) ••• A g n (x) 

Mo/o(«o) ••• Mo/m(^o) Xogo(uo) ••• Xog n (uo ) 


- hN foi^N ) ••• PNfm(uN) Ajvgo(wjv) Ajv5n(^iv)- 


( 2 . 6 ) 


By decomposing the top row as 


T (fo(x), ■ • • , f m (x), 0 , ■ ■ ■ , 0 ) + A(0, • • • , 0 , g 0 (x), ■■■ , g n (x)), 


(2.7) 


we have 

Pm,n (*^5 ^5 AO = Pm{P) ^ Qn{P)- 

On the other hand, the determinantal expression of Rm, n {x\ A, fi) implies 


Pm 1 n{y , k) ^k•> H>k) — H'k.P^kQni^k) — 0 — 0, 1, 2 , iV), 


( 2 . 8 ) 


(2.9) 
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and hence 


Pmi^k) • Qni^k) — ^k • l^k — 0 , 1 , . . . , N') 


( 2 . 10 ) 


as desired. 

The (IV+ 2) x (IV+ 2) determinants (12. 411 . (12.5)) representing P m (x) and Q n (x) can be condensed 
into an (m + 1) x (m + 1) and (n + 1) x (n + 1) determinants respectively, by means of a variation 
of the Dodgson condensation (see Appendix A). 

We denote by 


F (/?( M *))o <i<N, 0<j<m’ {Sj( U i))o<i<N, 0 <j<n (2-H) 

the matrices defined by the values of the functions fj(x) (0 < j < m) and (jj (x) (0 < j < n), 
respectively, at the reference points Ui (0 < i < N). We assume that the configuration of reference 
points Uk (k = 0,1,... , N) is generic in the sense that the minor determinants 

detF o ;.'.'.',n +n (0 <i<m), det(0 < i < n) (2.12) 


of maximal size with consecutive rows are all nonzero. 

Assuming that ^ 0, / 0 for k = 0,1,..., N, we set 


Ai 

Ah 


U id = —- det 


yf/jOi) 90 (Ui) 

1 ) go(u i+1 ) 


Aj+ n +i fj( u i+n+l) 9o( u i+n+l ) 
for 0 < i < m, 0 < j < m and 


Vij = ^ det 
^i 


9j{ui ) f 0 {ui) 

J^9j(u i+ 1 ) fo(u i+ i) 

/ij+TO+l 9j ( u i+m+ 1) /o( w i+m+l) 


5n(«i) 

Sn(^*+1) 

*?n(^i+n+l) 


:detG^;;;;^ +1 ) 1 


(2.13) 


fm{ u i ) 
/m(^i+1) 

fm(u i+m+1 ) 


: d etF 0 l + 1 :;^ +m+1 ) 1 (2.14) 


for 0 < i < n, 0 < j < n. Then by Lemma IA.2I (of condensation along a moving core), the 
(N + 2) x (N + 2) determinants (12.41) . (12.51) are condensed as follows into (m + 1) x (m + 1) and 
(n + 1) x (n + 1) determinants respectively (see also (IA.18I) . (IA.19I) L 


Theorem 2.2 The two functions P m (x), Q n (x) defined in Theorem \2.1\ are expressed as follows 
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in terms of (m + 1) x (m + 1) and (n + 1) x (n + 1) determinants respectively. 

AW --- fm(x) 

II II .”+». 


m— 1 n 

P m {x) = U Hi n A m+i det G o'X'.’X” det 

2=0 i =0 


Uqa 


Uq. 


n —1 m 


Qn(x ) = e m .n JJ Aj det F det 


2—0 2—0 


- 1,0 
go(x) 
Vo,o 

V n -\fi 


Um—l.m 
9n(x) ' 

V 0 ,n 

Un—l,n - 


(2.15) 


(2.16) 


where e m , n = (-1)™+™+". 

By expanding the determinants and along the first column we further obtain the series 
expansions 


n +1 


Uij = ^(-!) 


fc hi+k Aj 


k=0 

771+1 


r« = E*- 1 ) 


Aj+fc Pi 

fc Aj+fc //j 
Pi+k A i 


fj{ u i+k ) 


det G! 


ii,...,i+fc,...,i+n+l 
0,1....,n 


det G, 


2 + 1 ,.... 7 + 77 + 1 

0,1,...,n 


Sj iPi+k) 


det p i ’-’ i+k ’-t+m+i 


0,1 ,...,m 


? i+l,...,i+m+l 

,777 


(0 < i < m, 0 < j < m), (2.17) 


(0 < i < n, 0 < j <n). (2-18) 


k= o --- detF o,i,.’.:,( 

Hence the problem to determine P m (x) and Q n (x) is reduced to the computation of minor deter¬ 
minants of the matrices F = ( fj{ui)) i • and G = ( gj{ui)) i ■. We remark that these formulas for 
P m (x) and Q n (x) hold universally for any choice of the functions fj(x) and gj(x). 


In Sections 3 and 4, we show that these expansion formulas (12.171) . (12.181) in fact give rise to 
hypergeometric series of various types for appropriate choices of the functions fj(x), gj(x), the 
reference points Uk and the prescribed values vk = ^k/Pk- 


3 Hypergeometric series arising from determinants 


We explain below how the series expansions (12.171) . (12. 18[) can be used for generating hypergeometric 
series. In this section we use the notation of shifted factorials 


(a) n = a(a + 1) • • • (a + n — 1) = T(a + n)/T(a) (n = 0,1, 2,...). 
As a typical example, we consider the two sequences of rational functions 


(a + x)j 


(c + x)j 


9 ’ (x) = wPfs y=0il ’ 2 ’- ) 

with four complex parameters a, b , c, d, to form a pair (P m (x),Qn(x)) of rational functions 


Pm(x) — 'y ^ Pm,’ 


j =0 


0 + X)j 
(b + x)j 


Qn{x) — y ' Qn,j 


j =0 


(c + X)j 
(d + x)j ' 


(3.1) 


(3.2) 


(3.3) 
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Taking an arithmetic progression Uk = u + k (k = 0,1,2,..., JV; N = m + n) of points in C, we 
consider the Pade interpolation problem 


1 m (^ T k) _ A/- 

Qni,H T ^0 hk 


(k = 0,1,... ,1V). 


(3.4) 


Theorem 3.1 Consider the Pade interpolation problem (13. 3|) . (I3.4|) for the rational functions fj(x), 
gj(z ) in (13.211 and the reference points Uk = u + k (k = 0,1,..., IV; N = m + n). Then the solution 
(P m (x),Q n (x)) of Theorem, \2.1\ is explicitly given by 


where 


P m (x) = 


n;u * (d - c) 

EUo ( d + u m+i)n + 


m— 1 n 


j J ^ m-\-i det 


7=0 i=0 


Q(x)-e H ZcHh-a)* 
— c m,n j-rm 


n—1 m 


rii=0 (P + u n+i)m ~ =Q i=Q 


| ^ Aj det 


fo(x) 

U 0 ,o 

Um— 1,0 

9o{x) 

V 0 ,o 

14-1,0 


fm{x ) 
Uo,m 

Um—l,m - 


9n(x) 

V 0 , n 


u n - 


n— l,7i _ 


y^_l_ 

T _ (a + uf)j ( — n — 1 )k(d + 'Ui+n)fc (o + Ui+j)k(b + + & Ti+k At 

(6 + -u.,;)j •+( fc! (d + rtj) fc (a + rtj) fc (6 + rt i+ j) fc \ i+k m' 

, _ (c + Uj)j *(+ (~m - l)fc(fe +-»t +m )fc (c + rti +J )fc(d + ^)fc At + fc //i 
(d T uf)j k\ (b uf)k (c T rti)fc(d T i^i+j)fc gi+k A* 


(3.5) 


(3.6) 


(3.7) 

(3-8) 


As we remarked in the previous section, the functions LA- in Theorem 12.21 are expressed as 


71+1 




fc =0 


k Ti+k Ai 
Aj+fc Hi 


fjiP'i+k) 


det Gj 


<7,...,7+fc,..., 7+71+1 
0,1,...,77 

det G o!’i!. , ;;,n +n+1 


Since Uk = u + k (k = 0,1,... , IV), we have 

, , , _ (q + ^i+fc)j _ (a + Uj)j (a + tt»++(fr + 

J (b + u,i+k)j (6 + Ui)j (6 + Ui+j)k(a + Mt)fe 

The (n + 1) x (n + 1) minor determinants of the matrix 

G = ,(C + Mt)j 


(3-9) 


(3.10) 


\ i (3-11) 

. (a + Ut)j 7 0<i<N, 0<j<n 

can be computed by means of a special case of Krattenthaler’s determinant formula [5] (see Ap¬ 
pendix B). In fact by (1B.3I1 . we have 


det G 


i,i+ 1,... 

0,1,...,77 


, 7+77 


J / ^i J ...,i+fc J ...,i+n-|-l 
acr <-r 0)1 


nr=+(d-+ 

rii=o( d u i+l)n 

Y['iLi ll (d - c)i (-1 ) fc (-n - l) fc (d + tf»+n)fc 
n^=l (d + u i+l)n ( d 


(3.12) 
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Hence Uij is computed as 


j T _ (^ + Ui)j ^ “I - ^z+n)fc (& H“ ^i)k l^i-\-k 

Uil = W^), t __ 0 


(3.13) 


Jj k _q T (® “1“ Ui) k {b + u i+j)k A i+k hi 

The corresponding formula for Vij is obtained by exchanging the roles of ( m,n ), (a, c) and (6, d). 

If we choose the prescribed values appropriately, the series Uij and V t j in Theorem 13.11 give rise 
to generalized hypergeometric series 


r+l-fr 

Consider the case where 


) Oir 


ao, ol i, 

Pi, ■■■■ Pr 


E 

k =0 


^k (^1 )/c ' ' ' (^r)fc ( %\ 


Vk hk (h) k -■ ■ (t r ) k \wJ 


(«o)fc(ai)fc • • • (ov)fc k 

(1 )k{Pl)k ■ ■ ■ ( Pr)k 


(k = 0,l,2,...,N) 


with complex parameters si,..., s r and ti,..., t r . Since 

t^i+k^i _ (^1 H“ Ofe ’ ’ ’ (^r “1“ Ofc ^ 


^i+kt^i (^1 "h ^)k ’ ' ' (^r “1“ i)k ^ Z ' 


Uij and Vij are determined as 


TT _{a + Ui)j 

Uij — (u , ^ r+4-*r+3 


(6 + u 


ijj 


T /■ _ (c + «i)i 

Vij — / > | \ r+4U r +3 

(d + Ui)j 


—n — 1, d + Ui +n , a + Uj+j, 6 + itj, fi + *,..., t r + i w 
d + u i ,a j ru i ,b + u i+ j , s\+i,... ,s r + i ’2 

—m - 1, 6 + w i+n ,c + Ui +j ,d + Ui, si + i,... ,s r + «. £ 
6 + Uj,c + Ui,d + u i+ j,ti + i,... ,t r + i 1 w 


If we choose the prescribed values 


v k = — = 


A k (b + u) k (c + u) k (si)fe • • • (s r ) k (z\ k 


hk (a + u) k (d + u) k (ii)fc • • • (t r ) k \wJ 


- 


{k = 0,1, 2,, IV), 


then and Vu are slightly simplified as 


tt ( a + u Pj 

Uij — ,, . r+ 3 -T r +2 

(b + Ui)j 

_ {C + Ui)j 

Vij — / j , ^ r+3-*V+2 


(d + u. 


i>3 


—n — l,d + Ui+n, a + u i+ j,ti + i ,..., t r + i m w_ 
c + Ui,,b + u i+ j,si+i,...,s r + i ’ 2 

—m — 1,6 + Ui +n , c + Ui + j, si + i ,..., s r + i z 
CL + U{, d + Ui-j-j , tl + i,. . . , t r + Z ’ 10 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


(3.18) 


(3.19) 


As for the Pade interpolation problem for the rational functions fj(x) and gj(x) as in (13.21) . one 
can construct another type of determinant formula for P m (x) and Q n (x ) involving hyper geometric 
series. 
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Theorem 3.2 Consider the Pade interpolation problem (I3.3j) . (13.4p for the rational functions fj(x), 
gj(z) in (13.21) and the reference points Uk = u + k (k = 0,1,..., IV; N = m + n). Then the solution 
(P m (x),Q n (x)) of Theorem \2.1\ is expressed as 


N 


Pm(x) — K m ,nfc, d )11 Aj det 


i=0 


fo(x) 

$ 0,0 


L $771—1,0 


fm(x) 


$ 


0,m 


l,m _ 


N 


Qnipt) — det 


2=0 


ffo(x) 


ffn(x) 




0,0 


L $n—1,0 


9, 


0 ,n 


$71— l,n . 


where 


^ _ (a + u)j {~N)k(d + u + N — 1 — i)k (c + u) k (a + Uj) k (b + u) k g k 

ij ~TbTuj]hr 


(l)fe(c + u-i)k {d + u)k (b + Uj)k(a + u) k \ k ’ 


^ _ (c + it)j (— N) k (b + u + N — 1 — (a + ii)fc (c + Uj) k (d + u) k X k 

13 {d + u)j ||g (l)fc(a + M-i) fc (6 + ii) fc (d +«j)fc(c +ii) fc R.' 

We remark that if the prescribed values are given by 


_ _ (^l)fc * * ’ (Sr)fc / <2 

k Tk (h)k ■ ■ ■ (t r ) k 


(k = 0,1,2,..., N), 


then d>jj and give rise to generalized hypergeometric series 


r+5-^ r+4 


_ (o + u)j 
13 (b + u)j 
_ (c + u)j 

ij ~ (d + u)j r+5 r+4 
If we choose 


—N, d + u + N — 1 — i,c + u,a + u + j,b + u,ti,... ,t r w 
c + u — i,d + u,b + u + j,a + u, s\,..., s r ’ z 

—N, b + u + N — 1 — i,a + u,c + u + j,d + u,s\,... ,s r z 
a + u — i,b + u,d + u + j,c + u,t\,... ,t r ’ w 


Afc (b + u) k {c + u) k (si) k ■ ■ ■ (s r ) k ( z\k 


v k = — = 


Tk (« + u)k{d + u) k ( h) k ■ ■ ■ (t r ) k \w 


(k = 0,l,2,...,N), 


then and \Hj,- are simplified as 


_ (a + m)j 

(b + u)j 

_ (c + U)j 
13 (d + u) 3 


' r+3-fr+2 


r • 3 I> ■ 2 


—IV, d + u + N — 1 — i, a + u + j, t\,..., t r w 
c + u - i,b + u + j,si,... ,s r ’ z 

—N,b + u + N — l—i,c + u + j,si,...,s r _ z 
a + u — i,d + u + j,t\,... ,t r 'w 


(3.20) 


(3.21) 


(3.22) 

(3.23) 


(3.24) 


(3.25) 


(3.26) 


(3.27) 
























In order to obtain the expression of Theorem 13.21 we first rewrite (12.41) as 



fo(x) 

fm(x ) 

0 

0 

N 

JJAt det 

i=0 

E5 

0 I 0 

3.l'< 

• Y~fm{u 0 ) 

Ao 

go{uo) ■ 

• gn(u 0 ) 

fe; 

fell fe; 

_1 

• Y^fm(u N ) 

An 

go{u N ) ■ 

gni'U'N) 

N 

JjAj det 
i=0 

/(&) 0 

F G 





We construct an (N + 1) x (N + 1) invertible matrix L = such that ( LG)ij 

i+j<N, and define M = (Mjj)” by = (LG) m+i j. If we set = LF, we have 


'l 


"/(*) 


~f(x) 0 


'/(*) O' 




= 


= 

0 

L 


F G 


$ LG 


_ <f>" M_ 


Hence, by taking the determinants of the both sides we obtain 


P m (x) = A 0 • • • A tv det 


fo(x) ■ 

• fm{x) 


$00 

^0m 

detM 



detL ’ 

_$ro-l,0 • 

l,m_ 



which will give formula (13.201) with K mtn (c, d ) = det Mj det L. In view of 


9j{ u k ) 


(C + U k )j 

( d 4 - u k ) j 


(d + u) k (c + Uj) k ( c + u)j 
(c + u) k (d + Uj) k {d + u)j ’ 


we recall the Saalschiitz sum 


3 F 2 


—N, d + u + N — 1— i,c + u + j 
c + u — i,d + u + j 


{d - c) N (-i - j) N 
(c + u — i) N (d + u + j) N ’ 


namely, 


N 


E 


(~N)k{d + u + N — 1 — i)k(c + 
(1 )k{c + u-i) k (d + Uj) k 



(d — c)jy (—i - J)n 
(c + u-i) N (d + u + j) N ' 


With this observation, we define the matrix L = (Ljj)^._ 0 by 

^ _ (- N)j(d + u + N - 1 - i)j (c + u)j 
(l)j (c + u — i)j (d + u)j 


(0 < i,j < N). 


(3.28) 
= 0 for 

(3.29) 

(3.30) 

(3.31) 

(3.32) 

(3.33) 

(3.34) 
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The we have 


N 


(LG),, = £ 


(-N) k (d + u + N - 1 - i) k (c + Uj) k (c + u)j 


k =0 


(1 )j(c + u-i) k 


(d + Uj)k {d + u)j 


(d-c) N (-i-j) N ( c + u)j 

{c + u — i) N {d + u + j) N (d + u). 


(3.35) 


by the Saalschiitz sum. In particular (LG)ij = 0 (i + j < N). The determinant of the matrix M 
is computed as 


II 'll 

de.M = (-l)« , )n(LG) K ^=(-l)(”i , )n 

3=0 j =0 

Also, the entires of = LF are expressed as 


(d — c)n(—N)n 


(c + u)j 


(c + u - N + j) N (d + u + j)N (d + u)j 


N 


Qij = J2 L iky-Mu k ) 


k =0 
N 

E 

k =0 


Afc' 


(-N) k (d + u + N - 1 - i) k (c + u) k Hk (a + u k )j 


(1 )k(c + u-i) k 


(d + u)k A k (b + u k )j 


N 


(a + u)j y-r ( —N)k(d + u + N — 1 — i)k (c + u)k (a + Uj)k{b + u)k l^k 


( 6 + ^fe= o 


(l)fc(c + u-i) fc (d + u)fc (6 + Uj) fe (a + u)fe A fc ' 
The determinant of L can be computed again by Krattenthaler’s formula: 


det L = det 


(d -{- u -{- N — 1 — i)j ^ t- r ( —N)j (c + u)j 

i,j= 0 j = 0 + 


(c + u - i)j 

N 


( -i)( W 2 ~ 1 )tt (c zFFLL 

J- -L (r 4- n — 'iV 


l )j (~N)j(c + u)j 


3=0 


(c + U-j) N 


(d + u) 


(3.36) 


(3.37) 


(3.38) 


The constant factor in (13.20|) is determined as K rn n (c, d) = det Mj detL. 


4 Three types of very well-poised hypergeometric series 

In this section we consider three classes of hypergeometric series 

( 0 ) rational ... ordinary hyper geometric series 

( 1 ) trigonometric ... basic (or (/-)hypergeometric series 

( 2 ) elliptic ... elliptic hypergeometric series 

corresponding to the choice of a “fundamental” function [x]: 

(0) rational : [x] = e c ° x2+Cl x (0 = 0) 

(1) trigonometric : [x] = e c ° x2+Cl sin( 7 rx/cc) (O = Zcj) 

( 2 ) elliptic : [x] = e c ° x2+ci cr(x| 0 ) (O = Zwi © ZW 2 ) 
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where <r(x|fl) is the Weierstrass sigma function associated with the period lattice fl = Zuq © Zcu2- 
It is known that these classes of functions [x] are characterized by the so-called Riemann relation: 
For any x, a, (5 ,7 € C, 


[x + a] [x — a] [f3 + 7] [/? — 7] + [x + /3\ [x - /3] [7 + a] [7 - a] + [x + 7] [x - 7] [a + f3\ [a — /3\ = 0. (4.1) 
By the notation [x © y\ = \x + y\ \x — y\ of the product of two factors, this relation is expressed as 

[x © a] [/3 © 7] + [x © /3] [7 © a] + [x © 7] [a © /3] = 0. (4.2) 


In what follows, we fix a nonzero entire function [x] satisfying this functional equation. 
Fixing a generic constant <5, we define the 5-shifted factorials [x]*, by 

[x]k = [x\s,k = [x] [x + <5] ■ ■ ■ [x + (fc - 1)<5] (k = 0,1,2 ,...). 


(4.3) 


Then we define the very well-poised hyper geometric series r +*,V r +± [ao; «i ■ • • a r \z\ associated with 
[x] by 


r+5^r’+4 


a 0 ; a 1 , 


E 

k =0 


[ao © 2kd] 


[®o]fc [®l]fc ‘ ‘ ‘ [®r]fc 


[a 0 ] [^]fc[^ © o-o ai]fc [5 © ao a r ]fc 


z k . 


(4.4) 


In this paper we use this notation only for terminating series assuming that a* is of the form —nd 
(n = 0,1, 2,...) for some i. When z = lwe also write 


r+5^r+4 


a0; ai 


, , u, r 


E 

k =0 


[ao + 2 kS\ 


[ao]fc [ai]fc ■ ■ ■ [ar]fc 


[ao] [^]fc[^ © ®o ai]^ ■ ■ ■ [<5 © ao a r ]fc 


(4.5) 


In this notation, the celebrated Frenkel-Turaev sum is expressed as 


10^9 


ao; ai, a2, a,3,04, as 


[<5 + ao]jv[<5 + ao — ai — a 2 ]jv [(5 + ao — ai — a 3 ]jv [<5 + ao — 02 — 03 ] 


N 


[<5 + ao — ai]jv [<5 + ao — 02]^ [5 + ao — 03] at [<5 + ao — ai — 02 — 03] at ’ 
under the balancing condition ai + • • • + 05 = 2ao + <5 and the termination condition 05 = — N5 
(N = 0,1, 2,...). (See for example [3], |3]-) 

We remark that, in the rational case where [x] = x and <5 = 1, the r+ ^V r+ i series defined above 
is expressed in terms of a r+ 2 F r+ \ -series: 

ao? 2 ao © 1? ai, • • • a r 
2 ®o? , • • •, b r 

where 6* = 1 + ao — a* (i = 1 ,..., r). Also, in the trigonometric case where [x] = e cx / 2 — e ~ cx ^ 2 , 


r+5^r+4 


ao; ai, 


— r4-2 F 


r+2-fr+l 


(4.7) 


r+5^r+4 


a 0; ai, 


1 a?- 


E 

k =0 


l-q tp (t 0 \q)k {ti-,q)k ■■■ (t r -q) k 
1 - t 0 (q; q)k(qto/h;q)k ■ ■ • (qto/t r )k 


— r+3bF r+ 2 


to,h, ...,t r ,q,s 


(4.8) 
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in the notation of very well-poised g-hypergeometric series [3], where q = e cS , U = e CXi (i = 

r—l 

0,1,... , r) and s = (qto)^~z/t\ ■ ■ -t r . We discuss below a class of Pade interpolation problems 
that can be formulated in an unified manner in the three types of very well-poised hypergeometric 
series. 


Taking the two sequence of meromorphic functions 


fj(x) 

9j ( x ) 


[ a ± x]j 

[b ± x]j 

t c ± x \j 

[d zb x]j 


[a + x\j [a — x\j 
[b + x\j [6 - x]j ’ 
[c + x\j[ c-x\j 
[d + x]j[d - x\j 


(j = 0 , 1 , 2 ,...) 


(4.9) 


and the reference points u k = u + k5 (k = 0,1, 2,...), we consider the Pade interpolation problem 


Pm (^k) _ A k 

Qn(^k) pk 


(k = 0,1,... ,N) 


(4.10) 


for a pair of functions 


P m (x) =Pm,ofo(x) +Pm,lfl(x) H-b Pm ,m /mix'), 

Qn(x) = q n ,o 9 o(x) + q n ,i9i(x) H-b q n ,n 9 n(x) (4.11) 


where N = m + n. The prescribed values v k = A k /p k (k = 0,1,2 ... , N ) will be specified later. 


Theorem 4.1 Consider the Pade interpolation problem (14.101) . (14.111) for the functions fj(x), gj(x) 
in (TOD and the reference points Uk = u + k5 (k = 0,1,..., N\ N = m + n). Then the solution 
{P m (x),Q n (x)) of Theorem \2.1\ is explicitly given by 


r> f„\ _ n f -p-j- -pT \ 

Pm{x ) — C n (c, d) y- j.,, r ~ t ^ Pi ^ A m +i det 


TI/=0 ^ u m+l]r 


2=0 2=0 


n, (— r> („ TT \ TT,, 

Qn(x) — £m,nCm (o, 6) r j i i J. J. J. ^Ahi+i det 


foix) 

U 0 ,0 

. Um—1,0 

9o(x) 

V 0)0 


n;=0 i b ± u n+l]r 


2=0 2=0 


L Pn—1,0 


fm(x) 

Uo,m 

Um—l,m - 


9n{x) 



U n - 


n—l,n _ 


where 

n 

C n (c, d) = (-1)^) H\d-c\ k [c + d + (k- 1 )6} k , e m , n = (—l) mra+m + n 

k=1 


(4.12) 


(4.13) 


(4.14) 
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and 


jj _ [a ± Uj\j y-j [2 Uj + 2kS\ [2Uj\ k [-(n + 1)+ [uj-d + 5] k [uj + d + nd\ k 
lj [b ± Ui\j [2 Ui] [6\ k [2ui + (n + 2)6\ k [ui + d] k [ui - d + (1 - n)6]k 

[uj - a + 6\ k [uj + a + j5] k [ui + b\ k {uj - b+ (1 - j)5\ k fJ, i+k Xj 
[ui + a] k [[ui - a + (1 - j)6\ k [ui -b + 5] k [ui + b + jS\ k A i+k /jLi ’ 

y _ jc± + 2M] [2u+ [-(m + 1)+ + - b + + [itj + b + mS\ k 

lj [d±Ui}j^ o [2 Ui] [5\ k [2ui + (m + 2)5] k [ui + b\ k [ui - b + (1 - m)5\ k 

_ [-»j - c + + c + j+K + - d+ (1 -j)S\ k Xj+k^i 

[Ui + c]k[[ui - c + (1 - j)S]k[ui -d + 5\ k [ui + d + j5\ k m +k \i' 


(4.15) 


(4.16) 


As before we consider the expansion 


71+1 

Uij = £(-!)* 


k =0 


x+,...,z+/c,...,z+tz+1 

/^#7+1,... ,2+71+1 

U 0,l,...,7l 


fj (^Z+/c) 


/^z+/c ^z 
^z+fc /^z 


(4.17) 


of the determinant of (12.1311 . In this case we have 

, , . _ [a ± •u»+fc]j _ [a ± Uj]j [uj-a + 6} k [uj + a + jS\ k [uj + b\ k \uj - b + (1 - j)+ 

j Ul+k [b ± «;++ [b ± + + a\ k [[ui - a + (1 - j)5\ k [in - b + 5} k [ui + b + jd\ k ' 

The (n + 1) x (n + 1) minor determinants of the matrix 


G = ( 9j(ui )) 


0<i<N, 0 <j<n 


[c + Ui]j A 

[d ± Ui]j J o <i<N, 0 <j<n 


(4.19) 


can be computed by means of an elliptic extension of Krattenthaler’s formula (see Appendix B). 
In fact, by (jB.141) we have 


det G 


Z,Z+1 

0 , 1 ,... 


.. .,Z+7l 
,71 


/ oZ,...,Z+fc,...,Z+7l+l 

Ciet ^r 0 ,l,...,7Z 


C„(c, c?) 


C7 n (c, d) 


nr = 0 m ± «*++ 

nk [ 2tti + i + 

IL=C+ u i+l+l\n 


k [2uj + 2k5] [2uj]k[-{n + 1)<+ [ttj - d + + d + n+ 

[2«i] [<5] fe [2r/i + (n + 2)<+ [m + d\ k [ui - d+ (1 - n)5] k ' 


(4.20) 


Hence Uij is computed as 


Uij = 


[a ± Uj]j 

[b ± Uijj 


71+1 

E 


k =0 


[2 m + 2 kS] 

[2 ui] 


[2 Uj\ k [-(n + 1)+ [uj-d + 5} k [uj + d + nS] k 
[6\k[2ui + (n + 2 )5\k [Ui + d\ k [ui - d + (1 - n)5] k 


[uj - a + 5\ k [uj + a + j5\ k [ui + b\ k [uj - b+ (1 - j)5] k n i+k Xj 
[ui + a\ k [[ui - a + (1 - j)S] k [ui - b + 5} k [ui + b + j5\ k A i+fc /+ 


(4.21) 


The corresponding formula for V tJ is obtained by exchanging the roles of ( m,n ), (a, c) and (b,d). 
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Consider the case where the prescribed values are specified as 


Vk = — = 


A k _ ( z \ k j-T [u-e s + 5] 


n J 


Hk [ u + es j k 

Then we obtain very well-poised series 

[a ± Ui\j 


{k = 0,1 


UlJ [6 ± Ui\ 


■ r+nVr +11 2 iif, —(n + 1)5, Ui — d + 5, Ui + d + n5, 


Vit = f±^k r+nVr+u 


[d ± Ui] 


Ui — a + 5,Ui + a + j6, Ui + b,Ui - b + {1 - j)S, Ui + ei,..., u, + e r 
2ui \ — (m + 1)5, Ui — b + 5, Ui + 6 + n5, 




z J 


Ui- c+5, Ui+c+jd, Ui+ d,Ui— d+(l - j) 5 , Ui — ei+ 5 , ... ,Ui—e r +5 — . (4.24) 


w 


When 


A k _ ( z\ k [u - a + 5] k [u + b\ k [u + c] k [u -d + 5] fc [u -e s + S] k 


v k = — = \ ~ 


Hk \w/ [u + a\k[u - b + S] k [u - c + 6] k [u + d] k [u + e s ] k 
we obtain simpler very well-poised hyper geometric series 
[a ± m ]j 


Ua — r , , i r+lohr+9 


[b ± u 


*J J 


2 Ui] —(n + 1)5, Ui — c + 5,Ui + d + nd, 


Ui + a + j 5 , Ui-b + (l- j) 5 , Ui + ei,...,Ui + e T 


w 


z 


[c ± Ui]j 


^ [d ± u. 


r+10^r+9 


i\] 


2 Ui\ —(m + 1)5, Ui — a + 5 ,Ui + b + n 5 , 


Ui + c + j 5 , Ui — d + (1 — j) 8 , Ui — ei + 5 ,... ,Ui - e r + 8 


w J 


(4.22) 


(4.23) 


(4.25) 


(4.26) 


Another type of determinantal expression for P m (x) and Q n {x) is formulated as follows. We 
remark that this type of determinant formulas has also been discussed in HU- In what follows, we 
use the notation 




^05 ^1j ’ ' * i 


[clq + 2 k8\ [&o]/c [ai\ k ' ' * [&r]/c 

[ao] [5]fe[5 + ao — ai\ k • • • [5 + ao — a r ]k 


(4.27) 


Theorem 4.2 Consider the Pade interpolation problem (14.101) . (14.111) for the functions fj{x), gj(z ) 
in (IT9D and the reference points u k = u + k8 (k = 0,1,..., IV; N = m + n). Then the solution 
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(P m (x),Q n (x)) of Theorem \2.1\ is expressed as 



fo(x) ■ 

• fm(x) 

N 

Pm{%) — ^ det 

• O 

O 

d’o.m 

i= 0 

. ^m-1,0 • 



N 

Qn{x) — ^) l^i det 

i =0 


go(x) . 

• • 9n[X) 

*0,0 

■ ■ *0,n 

*n-1.0 • 

*n-l. 


(4.28) 


(4.29) 


where 


[a ± u]j 


N 


&ij = 7T"T— —NS, u — c + d +id, u + d + (N — 1)5 — id, u + c, u — d + 6, 
^ ± u w' k =o L 


u + a + jd, u — b + d — jd, u — a + d, u + b 


kk_ 

j Xk 


(4.30) 


= 


|- ky 2 u; —Nd, u — a + d + id,u + h + (N — 1)5 — id,u + a,u — b + 5, 

\d±u\j z —' 

L u fc=o 


AT 


ix + c + jd, u — d + 5 — j5, u — cd, u + d 
Consider the case where the prescribed values are specified as 


Afc 

J hk 


Vk = — = 


A k ( z\ k T1 [u - e s + d]k 


n 


Hk \wJ [u + e s ] k 


{k = 0,1,... ,1V). 


(4.31) 


(4.32) 


Then we obtain very well-poised series 


®ij = 


= 


[a ± u\j 
[b ± u\j 

[c ± u\j 
[d ± u] 


•+ 1414+13 


2 u; —Nd, u — c + d + id,u + d + (N — 1)5 — id, u + c, u — d + 5, 


u + a + jd,u — b + d — jd,u — a + d,u + b,u + e\,... ,u + e 1 


w 


z J 


"r+1414+13 


2 it; — 1V5, n — a-t-5 + z5, rt + 6-1- (IV — 1)5 — id, u + a, u — b + 5, 


u+c+jd, u—d+d—jd,u—c+d, u+d, u—e\ +5,..., u—e r +d 


w J 


(4.33) 


(4.34) 


When 


A k _ ( z\ k i u ~ a + 5]k[u + b\k ^ + c]k[u -d + 5]fc -pr [u - e s + 5] fc 
fJ-k Aw/ [u + a\ k [u-b + d] k [u-c + d} k [u + d\ k jj^ [u + e s \ k 


(4.35) 
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we obtain simpler very well-poised hyper geometric series 


_ [a±u]j 

Vi i — —-i r+10 lr+9 


^ lb±u lJ 

_ [c ± u]j 
[d ± u\j 


w 
z J 


2 it; — N5, u — c + d + id,u + d + (N — 1)<5 — id. 
u + a + jd, u — b + 5 — jd, u + e\ ,..., u + e r 
r _)_ioK -+9 2rt; — Nd, u — a + d + id,u + b + (N — 1)<5 — id, 

£ 

u+c+jd, u—d+d—jd, u—e± +5,..., u—e r +d 


wJ 


(4.36) 


(4.37) 


Theorem 14.21 can be proved by a procedure similar to the one we used in the previous section. 
In this case we define the matrix L = (Ljj )^._ 0 by 


Lij = 


2 u; —Nd, u — c + (1 + i)d, u + d + (.N — 1 — i)d, u + c,u — d + d 

for 0 < i,j < N. Then one can show 

[c + d + (j-i - l)d] N [-(i + j)d] N [2u + d] N [d - c\ N [c±u]j 


(• LG)ij = 


[u + d + jd]N[u - c + (1 - j)d]isr[u + c - id]N[—u + d - (1 + i)d ]n [d ± u]j 


(4.38) 


(4.39) 


by means of the Frenkel-Turaev sum, and hence ( LG\j = 0 for i + j < N . Then the series Ty are 
obtained by computing the product LF as before. We remark that in this case 


[c T d — (N + 1 — 2j)d]Ar[-IV<5]jv[2u + <5]jv[a! - c]n 


n 


j=o [ u+d+jd] N [u-c+(l-j)d] N [u+c-(N-j)d] N [-u+d-(N+l-j)d] N 


(4.40) 


The determinant of L can also be computed in a factorized form by the elliptic version (IB.14D of 
Krattenthaler’s formula: 


det L = 


n / =i[ ( %[ c + d+ (IV — 1 — 2 j)d]j[c-d — (N — 1 )d\j[2u +jd\j 


n£o[“ + c - i5\N[u -d-(N - 2 - z)<5]iv 

[2 u + 2 jd] [2 u\j [— Nd\j [u + c]j [u — d + d\j 

[2 it] [<5] j [2 u + ( N + 1)<5]j [u — c + d]j [u + d]j 


N 

n 

3= 0 


(4.41) 


The constant in (14.281) is given by K rn n (c\ d) = det M/ det L. 
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A Condensation of determinants 


In this Appendix A, we give a review on the variation of Dodgson condensation (Sylvester identity) 
of determinants due to Gasca-Lopez-Carmona-Ramirez [2], which we call the condensation along 
a moving core. For further generalizations of Sylvester’s identity, we refer the reader to Miihlbach- 
Gasca [9j. 

We first recall a standard version of the Dodgson condensation (Sylvester’s identity) for com¬ 
parison. For a general m x n matrix X = Kj -< n (with entries in a commutative 

ring), we denote by X^’'"’^ = (xi a j b )i <a<r l<b<s the r x s submatrix of X with row indices 
i \,..., i r £ {1,... , m} and column indices j i,..., j r £ {1,..., n}. When r = s, we denote by 
det A') 1 the corresponding minor determinant. 


Lemma A.l (Dodgson condensation, Sylvester’s identity) Let X = annxn ma¬ 

trix and set n = r+s (r, s > 1). We define anrxr matrix Y = {yij) ■ - =i by using the (s+1) x (s+1) 
minor determinants = det of X. Then the determinant ofY is expressed as 


det Y = det X (detX r r + 1 ;;;^)’- 1 ; Y = VH = X^+i;;;;;. 




(A.l) 


Proof : Define an n x n upper triangular matrix Z = (Zy)"' ._ 1 by setting 


z ij 


' Sij det X r r fi ^ 

< (-iy- r detX^+ 1 ; 




n 


(1 <i,j < r) 

(1 < i < r; r + 1 < j < n) 


■jij 


(otherwise). 


(A.2) 


Then for i = 1,..., r, the (i, ^-component of the product ZX is given by 


(ZX)ij - Zy { 


x 


ij T ^ ^ %ik %kj 

k=r -\-1 


n 

= detx;+ 1 ;;;; : >^+ £ {-i) k ~ r det x^ti’' 

k=r +1 


%kj 


= det X 




j,r+ 1, 


..,n 

..,n’ 


(A.3) 


namely, 


(^)y 


yij = 

0 (j = r + l,...,n) 


This means that 


ZX 


y o 

yr+l,...,n Tyr+l,...,n 


(A.4) 


(A.5) 
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Since det Z = zn ■ ■ ■ z rr , we obtain 


det X (det X;+l^Y = det Y det 

This implies the polynomial identity 

det Y = det X (det X;+ 1 1 ;;"’”) r - 1 

in the variables Xij (1 <i,j< n). 

When (r, s ) = (n — 1,1), (1A.1D means that 

det ( XijX nn x r!t x u j^j^ ^ | — det X x nn . 

Another extreme case (r, s) = (2, n — 2) implies 

det det - det det = det xtf;;;;* det X*;;;;£, 

which is equivalent to 

det X };;;;£ } det xi;;;;;2 - det 1 det xf;;;;;” , = det x$;;£ det X*;;;;;^. 


(A.6) 

(A.7) 

□ 

(A.8) 

(A.9) 

(A.10) 


These identities (1A.9I) . (1A.10I) are often referred to as Jacobi’s formula or Lewis-Carroll’s formula. 

The variant of Dodgson condensation that we use in this paper is the following identity due to 
Gasca-Lopez-Carmona-Ramirez |2]j. 


Lemma A.2 (Condensation along a moving core) Let X = (x^)(7_ 1 an n x n matrix and 
set n = r + s (r, s > 1). We define an r x r matrix Y = (yij) i - =1 by using the (s + 1) x (s + 1) 
minor determinants y = det of X. Then the determinant ofY is expressed as 


r —1 


det Y = det X det X\ 




r-WL, v‘i = ’• 


(All) 


i =1 


Proof: We define an n x n upper triangular matrix Z = x as follows by using s x s minor 

determinants of X: 


(-l)i 


~ l det vV-jV-A+s 
oer^ r+li n 




Then for i = 1,... ,r, we have 


(1 < i < r; i < j < i + s) 
(otherwise). 


i-\-s 


i-\-s 


( ZX ) = Y,z*x kj = ^(-I) fe -Metx, 


r+1 ,n 


x kj = del .v11 


k=i 


k=i 


(A.12) 


(A.13) 
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1 ... j 


r r + 1 


n 


Vij 


detX^ 


.,n 



Figure 1: Condensation along a moving core 


namely 


This means that 



(1 < j < r) 

(r + 1 < j <n). 


ZX 


Y 0 

yr+l,...,n yr+l,...,n 


(A.14) 


(A.15) 


Since det Z = z\ \ ■ ■ ■ z rr . we obtain 


det X J] det X l ;-\ l +* n = det Y det X r r +l’-’^. (A.16) 

i =1 

This implies the polynomial identity 


r —1 

det y = det A det X l r '^+* n 
2—1 


in the variables Xjj (1 <i,j <n). 

We remark that, if we renormalize the matrix Y by setting 


y = GW m =<“ x&fcT ( det 




(i,j = !,••• T), 


then equality (IA.11I) is rewritten equivalently as 


det A = det 1" det A£f};;;;£ 


(A.17) 

□ 


(A.18) 


(A.19) 
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B Variations of Krattenthaler’s determinant formula 

In this Appendix B, we recall Krattenthaler’s determinant formula j5] and its elliptic extension 
due to Warnaar [T5] , Although these formulas can be proved in various ways, we remark here 
that they are consequences of Lemma IB. 31 below, which can be regarded as an abstract form of 
Krattenthaler’s determinant formula (for recent works on the evaluation of determinants involving 
shifted factorials, see Normand (10]). 

We first recall a typical form of Krattenthaler’s determinant formula [5]. 

Lemma B.l For any set of variables Xi (0 < i < m) and parameters a k , /3k, 7fc, 5k (0 < k < m), 
one has 


det 


n 

0<k<j 


OtkXj + Pk 
'Yk^i &k 


m 


i,j =0 


][] (Xj - Xi) JJ {a k 5i - Pk'yi) 
0<i<j<m 0 <k<l<m 

n n (7 kXi + 5 k ) 

0 <i<m 0 <k<m 


(B.l) 


By specializing the parameters a k , /3 k , 7fc, 5 k , we obtain various determinant formulas. We quote 
below some of them. 

(a) Case where ak = 7& = 1 and Pk = a k , 5 k = b k : 


det 


n 

1 0 <k<j 


Xi + a k 

Xi + b k 


i,j =0 


] I (Xj - Xi) (bi - a k ) 

0<i<j<m 0 <k<l<m 

n n (*i+ 

0 < 2 <m 0 <k<m 


In particular, by setting a k = a + k, b k = b + k one obtains 

m 


/ (a + Xi)j \ m 

1(6 + 47^0 


(Xj - Xi) (6 - a)k 

0<i<j<m k =1 

J (P 3Ci)m 

0 <i<m 


where (a) k = a(a + 1) • • • (a + k — 1). 

(b) Case where /3 k = 5 k = 1, a k = -a k , 7fc = ~h- 

II _ 

Q<i<j<m 

n n 


det 


ni 

l 0 <k<j 


1 CLk'X'i 
- b k Xi 


I I {h - a k ) 

0 <k<l<m 
(1 - b k Xi) 


0 <i<m 0 <k<m 


(B.2) 


(B.3) 


(B.4) 
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By setting a k = p k a and b k = q k b, one has 


det 



0 axp, p)j 

( bxi;q)j 


m 


i,j =0 


(Xj-Xi) ( q l b-p k a) 

0 <i<j<m 0<k<l<m 

][ ( bx i - 1 q) m 

0 <i<m 


where ( a]p) k = (1 — a)(l — pa) ■ ■ ■ (1 — p k 1 a). In particular, 


det 



(' axf,q)j 

( bxi-q)j 


m 


i,j =o 


a(” +, )«(” + ‘) n te-x,)nwo;?), 

0 <i<j<m k= 1 

]Q ( bxi\q) m 

0 <i<m 


(c) Noting that 


set 


(az;q) n (ac/z;q) n = (1 + a 2 cq 2k - aq k (z + c/z)) 

0 <k<n 


(B.5) 


(B. 6 ) 


(B.7) 


Xi = Zi+c/zi; a k = -ap k , /3 k = 1 + a 2 cp 2k , - y k =-bq k , 5 k = l + b 2 cq 2k . (B. 8 ) 


Then we have 

det ( (az i -p) j (ac/z i ;p) j \ m 

\{bz i ;q) j (bc/z i -,q) j ) ij=0 

(Zj - Zi)(l - c/ziZj) {bq l - ap k )(l - p k q l abc) 

0 <i<j<m 0 <k<l<m 

II ( bz i ;«) m ( bc/zi-q ) m 

0 <i<m 

In particular, 

det / {azi',q)j{ac/z i ;q) j \ m 
V (t>z i ;q) j {bc/z i ;q) j ) iJ=Q 

m 

t) q { 3) ( Zi - Zj)(l - c/ziZj) Yl(b/a;q) k (q 2 ( - m ~ 1 - k '>abc;q) k 

0 <i<j<m k =1 

n m ( bc/zi-q ) m 

0 <i<m 


(B.9) 


(B.10) 


Let [x] a nonzero entire function in x € C and suppose that [x] satisfies the so-called Riemann 
relation: For any x, a,f3 ,7 € C, 

[x ± a] [f3 ± 7 ] + [x ± /3 ][7 ± cc] + [x ± 7 ][a ± /?] = 0, (B.ll) 
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where [x ± a] = [x + a] [x — a]. This functional equation is equivalent to 

[x ± u] [y ± v\ — [x ± v\ [y ± u] = [x ± y] [u ± v]. 


(B.12) 


The following lemma is the elliptic extension of Lemma IB. II due to Warnaar [15]. 

Lemma B.2 For any set of variables Xj (0 < i < m) and parameters ak,bk (0 < k < m), one has 


det 


n 

0<k<j 


[Qfc ± Xi] 
[b k ± Xi] 


m 


hj =o 


n ± ^ n ^ 

0 <i<j<m 0 <k<l<m 

n n ± Xi i 

0<z<m 0 <k<m 


As a special case where ak = a + k5, bk = b + kd (0 < k < m), we obtain 


/ M4 \ ro 

\[b± Xi\j ) ij=0 


JJ [x i ±x j ]Y]_[b~ a] k [a + b +(k - l)6] k 

0<i<j<m k =1 

| \p i ^z]m 

0<i<m 


(B.13) 


(B.14) 


where [a]k = [a][o + <5] • • • [a + (fc — 1 )<5] and [a ± b\k = [a + b\k[a — 6]^. 

Lemma IB. II and Lemma IB. 21 can be proved as consequences of the following abstract form of 
Krattenthaler’s determinant formula. 

Let aik, bik (0 < i < N] 0 < k < N ) be elements a field IK with bik ^ 0 for all i, k, and consider 
the matrix 


X m = 



(B.15) 


for m = 0,1,..., N. Suppose that there exist elements Pij (0 < i,j < N), qki (0 < k < l < N ) of 
IK such that 


O'ikhjl O’jkbil — PijQkli Pji — Pij 
for all*, j G {0,1..., N} and k, l € {0,1,..., N — 1}. 

Lemma B.3 Under the assumption (IB.16I) . the determinant det X m is factorized as 


det X m 



11 Pji II qkl 

0<i<j<m 0 <k<l<m 

m m —1 

n n b hk 

i =0 fc=0 


for m = 0,1 ,... , N. 


(B.16) 


(B.17) 
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Set T m = det X m for m = 0,1, 2,... , N, so that 


to 


h 


ti = det 


1 

1 


qqo ~ 

^00 

QlO 

bio- 


Oio _ «00 
&10 &00 


The hrst nontrivial case is guaranteed by the assumption (IB.161) : 


Tl 


Qio^oo ~ apo^io 
boobio 


P 10 goo 
boo b io 


(B.18) 


(B.19) 


Then (1B.17|) can be proved by means of the Lewis-Carroll formula. In fact from (IA.10I) . we obtain 
the bilinear identities 


^TTL- (— 1,1 rri rri f \ 

T T m J-C-L R\Tm) 

0771+1,1 


an 

bn 


Tc{Tm)T R (T m ) 


T m +1 TcTR{T m - 1), 


(B.20) 


for r m , where Tr and Tq stand for the symbolic shift operator for the row indices and the column 
indices: 


Tr( p%j ) — T+1 .j: (bjj j — b l+ i j , T]i(pij ) —i?i+i,j+i; T]i(qki) — 

Tcfaij) = Ojj+i, Tc(bij) = hj + i, Tc{pij)=Pij, Tc(q k i) = 9fc+i,H-i- 

Thanks to the bilinear identities, one can prove 


(B.21) 


/ 7~m — 


n ph n qki 

0<i<j<m 0 <k<l<m 


m m— 1 

nir 


(m = 0,1,... ,1V) 


(B.22) 


i.k 


i =0 k =0 

by the induction on m. 

We remark that Lemma [Bj] is the case where a,+ = + fa, bik = i k Xi + S k . Since 


+ fa)(llXj + 5i) - (afcXj + Pk)(liXi + Si) = (xj - Xj)(a k 6i - fall) (B.23) 

the factorization condition (IB.161) is verified with p VJ = Xi — Xj and q k i = a k Si — fall- Lemma IB. 21 
is the case where = [a k ± Xj \, b, k = [b k ± Xj \. Since 

[a k ± Xi] [bi ± xj] - [a k ± xj] [bj, ± x;] = [a k ± b{\ [x t ± xj] (B.24) 

the condition (IB.16D is satisfied with pij = [x* ± Xj] and q k i = [a k ± bj\. One can prove in fact that 
generic solutions to the system of equations (IB. 161) reduce to the case of (IB.23D . It would be worth¬ 
while, however, to recognize the role of bilinear equations (IB.211) which lead to the factorization of 
determinants. 


23 




























References 


[1] A.G. Akritas, E.K. Akritas and G.I. Malaschonok: Various proofs of Sylvester’s (determinant) 
identity, Math. Comput. Simulation 42 (1996), 585-593. 

[2] M. Gasca, A. Lopez-Carmona and V. Ramirez: A generalization of Sylvester’s identity on 
determinants and the application to interpolation problems, in Multivariate approximation 
theory, II (W. Schempp and K. Zeller, Eds.), pp. 171-184, Birkhauser, 1982. 

[3] G. Gasper and M. Rahman: Basic hypergeometric series , with a foreword by Richard Askey, 
Second edition. Encyclopedia of Mathematics and its Applications, 96 , Cambridge University 
Press, 2004. xxvi+428 pp. 

[4] Y. Kajihara and M. Noumi: Multiple elliptic hypergeometric series: An approach from the 
Cauchy determinant, Indag. Math. 14(2003), 395-421. 

[5] C. Krattenthaler: Advanced determinant calculus, The Andrews Festschrift (Maratea, 1998), 
Sem. Lothar. Combin. 42 (1999), Art. B42q, 67 pp. (electronic). 

[6] C. Krattenthaler: Advanced determinant calculus: a complement, Linear Algebra Appl. 411 
(2005), 68-166. 

[7] Y.L. Luke: Mathematical functions and their approximations , Academic Press, Inc. [Harcourt 
Brace Jovanovich, Publishers], 1975. 

[8] Y.L. Luke: Algorithm for the computation of mathematical functions , Academic Press Inc. 
[Harcourt Brace Jovanovich, Publishers], 1977. 

[9] G. Miihlbach and M. Gasca: A generalization of Sylvester’s identity on determinants and some 
applications, Linear Algebra Appl. 66(1985), 221-234. 

[10] J. M. Normand: Calculation of some determinants using the s-shifted factorial, J. Phys. A 37 
(2004), 5737-5762. 

[11] M. Noumi, S. Tsujimoto and Y. Yarnada: Pade interpolation for elliptic Painleve equation, 
Symmetries, integrable systems and representations , 463-482, Springer Proc. Math. Stat., 40 , 
Springer, 2013. 

[12] H. Rosengren and M. Schlosser: Summations and transformations for multiple basic and elliptic 
hypergeometric series by determinant evaluation, Indag. Math. (N.S.) 14 (2003), 483-513. 

[13] P.E. Spicer, F.W. Nijhoff and P.H. van der Kamp: Higher analogues of the discrete-time Toda 
equation and the quotient-difference algorithm, Nonlinearity 24(2011), 2229-2263. 


24 



[14] V.P. Spiridonov and A.S. Zhedanov: Elliptic grids, rational functions and the Pade interpola¬ 
tion, Ramanujan J. 13(2007), 285-310. 

[15] S.O. Warnaar, Summation and transformation formulas for elliptic hypergeometric series, Con- 
str. Approx. 18 (2002), no. 4, 479-502. 

[16] Y. Yarnada: Pade method to Painleve equations, Funkcial. Ekvac. 52(2009), 83-92. 

[17] Y. Yarnada: A Lax formalism for the elliptic difference Painleve equation, SIGMA 5(2009), 
042 (15 pages). 


25 



